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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 


TECHNICAL NOTE NO. 1876 

CALCULATION OF THE AERODYNAMIC LOADING OF FLEXIBLE WINGS 
OF ARBITRARY PLAN FORM AND STIFFNESS 
By Franklin W. Diederich 

SUMMARY 


A method is presented for calculating the aerodynamic loading, the 
divergence speed, and certain stability derivatives of wings and tail 
surfaces of arbitrary plan form and stiffness. Provision is made for 
using either stiffness curves and root— rotation constants or influence 
coefficients in the analysis. Computing forms, tables of numerical 
constants required in the analysis, and an illustrative example are 
included to facilitate calculations by means of the method. 


INTRODUCTION 


The distribution of the aerodynamic loading on wings and tail 
surfaces is important both for the structural analysis of these com- 
ponents, since it determines the applied bending moment and torque 
acting at any station, and for their aerodynamic analysis, since it 
affects the stability derivatives to a large extent. At high speeds the 
aerodynamic loading, particularly in the case of swept wings, is 
greatly affected by the structural deformations caused by the loading. 

The present paper is concerned with the determination of the effects of 
structural flexibility on the aerodynamic loading of wings of arbitrary 
plan form and stiffness. 

The problem of load distribution was analyzed for unswept flexible 
wings as early as 1926 (reference l) but has received relatively little 
attention since that time. The only new effect considered in subsequent 
work is -aerodynamic induction (reference 2). No work appears to have been 
done on the loading of flexible swept wings. The related problem of aero- 
elastic divergence of swept wings with certain prescribed stiffness 
variations has been treated in reference 3* 

The present paper treats the problem of aerodynamic loading by 
matrix methods. Aerodynamic induction is taken into account approxi- 
mately, since suitable aerodynamic influence coefficients are not avail- 
able for wings of arbitrary plan form. When they become available they 
can readily be incorporated in this method. Structural flexibility is 
taken into account in the form of either calculated stiffness variations 
or measured influence coefficients. The required integrating matrices 
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are presented for "both a six-point and a ten— point solution. For the . 
six-point solution convenient computing forms are included as well. The 
method is illustrated by means of an example . In addition to the analysis 
of the aerodynamic loading, the determination of the related divergence 
speed and of certain stability derivatives is discussed. 

For the convenience of the reader unfamiliar with matrix terminology 
a summary of matrix methods has been included in the appendix. The 
sections on " Application of the Method" and, in particular, "Instructions- 
for Solution" may be re^id without reference to the section "Derivation of 
the Method." 


SYMBOLS 


A 

[a] 


a* 

ac 

b 

c 






V 


‘H 


- © 


aspect ratio 


aeroelastic matrix 


dimensionless parameter 


n e^ s A^ e l r c r^ cos ^ A 

(GJ) 


parameter (m e qs w c r cos a) 

section aerodynamic center, measured from leading edge, 
fraction of chord 

wing span, inches 

chord measured parallel to the air stream, inches 
average wing chord, inches 

section lift coefficient (—-] 

■ \qc 


wing lift coefficient 




wing root bending-moment coefficient 


'M» N 
s qSb, 


\ 


wing rolling-moment coefficient (bCgj^) 
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El bending stiffness in planes perpendicular to the elastic 

axis, pound— inches 2 

e location of elastic axis measured from leading edge, fraction 

of chord 



distance from reference axis to section aerodynamic center 
(positive forward), fraction of chord 

torsional stiffness in planes perpendicular to the elastic 
axis, pound— inches^ 

unit matrix 

matrix defined by equation ( 12 ) 

integrating matrix for single integration from tip to root 
first row of matrix 

integrating matrix for double integration from tip to root 
first row of K2 matrix 

integrating matrix for single integration from root to tip 

matrix relating concentrated and accumulated torque 

matrix relating concentrated loads and accumulated bending 
moments 

matrix converting torques due to distributed loads to torques 
due to concentrated torques 

matrix converting bending moments due to distributed loads to 
bending moments due to concentrated loads 

lift on both wings but excluding lift on part of wing covered 
by fuselage, pounds 

running air load along the reference axis, pounds per inch 

accumulated bending moment (in planes perpendicular to the 
reference axis unless specified otherwise), inch— pounds 

effective section lift-curve slope for angles of attack due to 
deformation, per radian 
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effective section, lift— curve slope for additional— type angle— 
of -attack distributions, per radian 

section lift-curve slope, per radian 

concentrated load, pounds 

wing— tip helix angle 


root— twist constants (see equation (9)) 
root-bending constant (see equation (9)) 
dy nami c pressure, pounds per square inch 
concentrated torque, inch— pounds 

total wing area including part of wing covered by fuselage, 
square inches 

distance from wing root to wing tip perpendicular to the air 
stream (see fig. l), inches 

length of wing along reference axis (see fig. l), inches 

accumulated torque (in planes perpendicular to the reference 
axis unless specified otherwise), inch— pounds 

distance between the effective root and the innermost complete 
section of the torsion box perpendicular to the elastic 
axis, inches 

lateral ordinate measured from wing root, inches 
lateral center of pressure, inches 


angle of attack, radians 
equivalent angle of attack, radians 



local dihedral angle due to deformation or slope of wing 
deflection curve at reference axis, radians 

structural deflection, inches 

distance along reference axis, inches 

angle of sweepback (measured to the referenoe axis unless 
specified otherwise), degrees 
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influence-coefficient matrix for wing twist in planes parallel 
to the air stream due to concentrated unit loads applied at 
the reference axis, radians per pound 

influence-coefficient matrix for wing twist in planes parallel 
to the air stream due to concentrated unit torques applied 
in planes parallel to the air stream, radians per inch- 
pound 

angle of twist in planes perpendicular to the reference axis, 
radians 


Subscripts: 


a 

c 

D 

fw 

g 

LE 

M 

MAC 

P 

r 

rw 

s 

sub 

spr 

T 

TE 

w 


additional 

midchord 

divergence 

flexible wing 

geometric 

leading edge 

due to bending moment 

pertaining to the mean aerodynamic chord 

damping in roll 

at root or effective root 

rigid wing 

structural (due to structural deformations) 

subsonic 

supersonic 

due to torque 

trailing edge 

wing exclusive of fuselage 


Prime mark: 


in or pertaining to sections parallel to air stream rather 
than perpendicular to the reference axis 


Matrix notation: 

{ | column matrix 

LJ 


row matrix 
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[] 

square matrix 

[°] 

diagonal matrix 


DERIVATION OF TEE METHOD 
Method Employing Stiffness Curves 


Assumptions. — In the development of the method the following 
assumptions are made: 

(a) The effects of aerodynamic induction may he taken into account 
by applying an over-all correction to strip theory and rounding off the 
resulting load distribution at the tip. 

(b ) All deflections and angles of attack are smal 1 . 

(c) The wing is mounted flexibly at an effective root perpendicular 
to the elastic axis through the intersection of the elastic axis an d 
the fuselage (see fig. l), the root rotations being proportional to the 
root bending moment and root torque. 

(d) An elastic axis exists in the outer portion of the wing, this 
axis being defined as the elastic axis the wing would have if it were 
mounted rigidly some distance outboard of the root approximately perpen- 
dicular to the midchord line. (Near the root the elastic axis is defined 
as the extension of the outboard elastic axis.) 

(e) All deformat ions are given by the elementary theories of 
bending and of torsion about the reference axis, which in this case is 
the elastic axis. 

Air loads .— In keeping with assumptions (a) and (b) the force on 
a wing section of unit width parallel to the direction of flight is 


2 * = 



+ J cos A 


= mQqca cos A 


( 1 ) 


where the equivalent angle of attack a is defined by 
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The effective section lift-curve Blope nig for angle-of -attack distribu- 
tions due to structural deformations has been given in reference 3 on 
the basis of the reasoning of reference k for subsonic speeds as 


mg = mo 


A 

A + k cos A 


( 2 ) 


4 

Similarly, the effective lift-curve slope m^ for additional-type 
angle— of -attack distributions is determined by the same reasoning as 



= “o 


A 

A + 2 cos A 


( 3 ) 


Both slopes must be multiplied by cos A, as in equation (l), in order 
to apply to loads acting on sections and angles- of attack measured 
in planes parallel to the direction of flight. 

The torque of the running load Z* about the reference axis is Z*e 1 c 
for uncambered sections (for cambered sections the torque at zero lift 
must be added and the analysis of the following paragraphs modified 
accordingly). This torque may be resolved into a running torque about 
the elastic axis and a running bending moment about a line perpendicular 
to the elastic axis. The running load, torque, and moment must then be 
multiplied by. cos A to yield their values per unit length along the 
elastic axis, so that 

Z = mgqca cos 2 A (lj.) 

or, in matrix notation, 

{l} = m e q(cos 2 A) [c]{a} (ka) 


, The running torque and the running bending moment are, respectively, 

jze^c cos A] and j— Ze^ sin A j. The running bending moment leads to 

accumulated bending moments which have to be added to the accumulated 
bending moment due to the running load. 

The accumulated torque T is obtained from the running torque by 
an integration inboard from the tip. This integration may be performed 
by a matrix JeJ which is based on Simpson* s rule with a modification 

suggested by V. M. Falkner at the tip. (See appendix. ) The effect of 
Falkner*s modification is to round off the calculated load distribution 
and cause it to go to zero with an infinite slope at the tip, as the 
aerodynamic lift distributions actually do. The matrix is given’ in 
table I. 
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Similarly, the accumulated tending moment M is obtained by a double 
integration inboard from the tip of the running load and a single 
integration of the running moment. The double integration may be 
performed by another matrix [Kg] (given in table II ), which is based 

on the equivalent of Simpson’s rule for moments, Falkner’s modification 
again being made at the tip. The derivation of the integrating matrices 
is .discussed in somewhat greater detail in the appendix. 

The accumulated torque and bending moment may then be written as 


M - 8 a[ k i] ( le •j^C CO8 A} 


and 


= m eT s A e l r c r 2 ( cos 3 a )[ k i] 




■a 2 ^] (*} - s a[%] 


{ 


ie^c sin 





( 5 ) 


(6) 


Equations of equilibrium .— The equations of equilibrium of a 
deformed wing referred to the elastic axis are 

( 7 ) 


( 8 ) 


These equations must be integrated outboard from the root to obtain <p 
and T. The integrations may be performed by a matrix (see 

table III and appendix), also based on Simpson’s rule without the tip 
modification, however, since the torques and moments go to zero with 
finite and zero slopes, respectively. To the deformations obtained in 
this manner the rotations due to the root deflection, cp r and f r . 


GJ | 2 = T 
dy 


dr 

El — = M 
dy 
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must be added. The root rotations are defined by four dimensionless 
constants: 


Vp 


W 


w/(GJ) 


(9a) 


w/(GJ) r 


(9b) 


Q r m = 


_ Fr r /' 


T 


W 


/(El), 


(9c) 


ft r 

M w/(El) 




(9d) 


which may be combined into two other constants 




%1 


orr^r 

L (Gj) r 

w/(OJ) r 

(> (El) r 

“nA 

/. <">r 

w/(Gj) r 

V* (El) r 


tan AQp^ I cos A 


*Qp. COS . 


'My 


(9e) 


(9f) 


w being defined as in figure 1. The deformations may then be written as 
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and 


( r } = 


“A 


(El) 


N 


(El), 


EE 


it 


( M } + ^«r T [io]W 


where the matrix is defined by 


( 11 ) 


w- 


0 0 0 0 
10 0 0 
10 0 0 
10 0 0 
10 0 0 


J 


(12) 


TJie angle of attack due to the structural deformations a 8 is 
related to <p and T "by 

a a = (9 — f tan A) cos A (13) 


If equations (5), (6), (l0), and (ll) are substituted in the matrix 
equivalent of equation (13), the following relation is obtained.' 


{a s } = a[A]|a| 

where the aeroelastic matrix [a] is defined by 


M- 


N 


O 

(GJ), 


GJ 


By (VP ^ tai1 ^ [l0] 




N 


fl /c 

e l. 


© 


m 




e-j c r cos A Sdk 


.] H [4] 


(15) 
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and the parameters 


a 


m e qsA 2 e 1 ^c r 2 cos4A 


(GJ) r 


(16) 


and 


d 

a 


My % 

(El ) r e^c r cos A 


tan A 


(l6a) 


are subst antially the same parameters as those used in reference 3* 

Solution of the equations .- If it is desired to calculate the 
aerodynamic loading corresponding to a given geometrical angle-of-attack 
distribution and dynamic pressure, equation (l4) may be rewritten as 
follows: 



( 17 ) 


In this form it constitutes a set of linear simultaneous equations for 
the ct, values in terms of a,g values. Prom the calculated cc values 

the lift distribution may be determined from either equation (l) or (4). 


. The divergence dynamic pressure may be obtained from equation (17) 
by setting the determinant of the square matrix on the left side of the 
equation equal to zero. This procedure is equivalent to setting a g 

equal to zero in the term a of equation (l4), so that 


{ a s} = a [ A ]{as} (18) 

The critical value of a is then determined by matrix iteration and 
hence the divergence dynamic pressure from equation (l6). 


Method Employing Influence Coefficients 

The assumptions made in the preceding sections concerning the 
behavior of the wing structure are unnecessary if influence coefficients 
for the given structure are available from test data or refined methods 
of calculation. The coefficients most convenient for this analysis are 
those giving the rotation of the structure in planes parallel to the 
direction of flight due to vertical loads applied along a convenient 
reference axis and due to torques about lines perpendicular to the 
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direction of flight. Since it is usually more convenient to apply 
concentrated rather than distributed loads in structural tests, the 
influence coefficients are considered in this analysis to have 
been obtained in this manner. 

The angle of structural deformation a g may be expressed in terms 
of the influence coefficients and $p as follows: 

{<%} - MW + MW (19) 


where the R*s and P*s are arbitrary concentrated torques and loads, 
the latter being applied at the reference axis. The accumulated torques 
and bending moments about lines perpendicular and parallel, respectively, 
to the direction of flight may be related to the concentrated torques 
and loads by means of the summation matrices and 

appendix) as follows: 


(see 


|t*} = [k 4 ]|r} - tan a{m»} 

M - svMW 


( 20 ) 

( 21 ) 


These relations may be solved for the values of R and P required to 
produce given distributions of accumulated torque and bending moment 


M= M + tan Ajl 

w ■ kbfh 


(22) 


(23) 


The accumulated torques and bending moments produced by the air 
load are then 


{t»} = s w[ K ij{* ,e l c } - |m«J tan A 


( 24 ) 


{M*} = s v 2 [K2]{l*} 


( 25 ) 
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Upon substituting equations (22), (23), (24), (25), and. (l) into 
equation (19), the following equation is obtained: 

(a 8 }= a*[A*][a] (2 6) 

where 


a* = mgqSyOj. cos A 

and 

M - %^NN ^ 


where, in turn 

[*7] * fs] K) 

are given in tables IV and V. 

The solution of equation (26) is obtained in the manner previously 
described for equation (14)., 


APPLICATION OF THE METHOD 
Dete r m in ation of the Structural Parameters 


At the time an aeroelastic a n alysis is performed no experimental 
stiffness data are usually available, so that the calculated stiffness 
curves must be used. In order to use these curves it is necessary to 
assume the existence of a reasonably straight elastic axis. The location 
of this axis may be estimated by considering it to be the line connect ing 
the shear centers of the individual sections. If the elastic axis 
obtained in this manner is not reasonably straight within a few percent 
of the chord, the results of the analysis may not be sufficiently reliable. 


(29a) 

(29b) 


(27) 

(28) 
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The stiffnesses GJ and .El do not have much physical significance 
inboard of the last point where there is a complete cross section of the 
torsion box. (See fig. 1.) In order to arrive at estimates of the root 
stiffnesses (Gj) r and (Ej) r , which serve primarily as reference values 

in this analysis, the stiffness curves have to be extended. It is con- 
venient to consider the stiffnesses to be constant inboard of the last 
complete section of the torsion box; this procedure should yield conserva- 
tive values of the root rotations. 

The most difficult problem incurred in analyzing the deflections on 
the basis of stiffness curves appears to be the estimation of the root 
rotations. As used in this analysis, they are the torsion and bending 
deflections imposed by the triangular inner portion of the wing and the 
carry-through bay on the rest of the wing. As seen in figure 2, which 
is plotted from the data of reference 5, these values are essentially 
constant along the span, so that they actually constitute rigid— body 
rotations. (The bending rotations have been obtained by talcing the 
difference in slope between curves calculated by considering the wing 
to be cantilevered at the effective root — the root used to calculate 
torsional deformations in reference 5 - and the averages of the leading- 
edge and trailing-edge deflections actually measured. The twists were 
obtained by subtracting the twists calculated on the basis of the assumed 
effective root from the measured twists.) 

The rotations should in any practical case be calculated by analyzing 
the triangular root and the carry— through bay and made dimensionless by 
means of equations (9)» If such an analysis is not available, the dimension- 
less rotation parameters shown in. figure 2 may be used as a guide; it 
must be kept in mind, however, that in the case of a sweptforward wing 
the parameters and Qp^ would have the opposite sign. 

Once the structure under investigation is built, fairly Bimple 
deflection tests, similar to those performed in reference 5, may be used 
to check the root— rotation parameters by calculating the differences 
between the observed rotations and those calculated by simple beam 
theory considering the wing cantilevered at the effective root; at the 
same time the existence and est ima ted location of the elastic axis may 
be verified. If the experl mental program is fairly extensive it is 
desirable to measure influence coefficients directly. These influence 
coefficients can then be used in conjunction with the alternate method 
described in the preceding section to obtain a quick check on the aero- 
elastic analysis based on calculated stiffnesses. 

The influence coefficients used in the analysis consist of the 
rotations of sections parallel to the direction of flight due to 
concentrated unit torques in planes parallel to the plane of symmetry 
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or concentrated unit loads at the reference line. These rotations in 
radians are entered in tables of the form: 


TWIST AT STATION 7/ s w DUE TO 
UNIT CONCENTRATED TORQUE AT 


KJ 

TWIST AT STATION j/b^ DUE TO 
UNIT CONCENTRATED LOAD AT 




yi/sw 

7/*v 

0.2 

0.4 

0.6 

0.8 

0.9 

1.0 

0 






■ 

0.2 






■ 

0.4 







0.6 







0.8 







0.9 








y 1K 


\ 

y/ew n. 

0.2 

0.4 

0.6 

0.8 

0.9 

1.0 

0 

■ 

■ 

■ 

■ 


■ 

0.2 



■ 

■ 

■ 

■ 

0.4 







0.6 







0.8 


- 





0.9 








These particular tables would be used for a six-point analysis; similar 
tables would be used for a ten-point analysis. In either case it is 
to be noted that the twists are measured at values of y / By from 0 to 0.9, 

whereas the loads are applied at y^/sy. values from 0.2 to 1.0. The 

tables obtained in this manner constitute the desired influence- 
coefficient matrices. 

If the wing sections are found to twist nonuniformly, so that they 
•become cambered in effect, the angles of twist a s to be entered in 
the influence-coefficient matrices have to be defined in a different 
mann er according to whether the aeroelastic analysis is performed for 
subsonic or supersonic speeds. At subsonic speeds the lift depends on 
the slope of the mean camber line at the three — quarter — chord point, so 
that the effective angle of attack is 
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At supersonic speeds the lift depends primarily on the average slope of 
the mean camber line, so that 


a B = 


^LE ^TE 


(31) 


Determination of the Aerodynamic Parameters 

The selection of the aerodynamic parameters and e 1 for the 

calculation of the divergence speed has been discussed in reference 3. 

For calculating the aerodynamic loading at a given flight condition the 
aerodynamic parameters are chosen for that flight condition. The 
effective lift-curve slopes i% and m ei are applicable only to subsonic 

subcritical speeds. At higher speeds no simple span correction is 
available; neglect of the span correction tends to be conservative for 
calculation of the divergence speed and the aerodynamic loading, however. 


Instructions for Solution 

Two sets of integrating matrices have been prepared, one for a six- 
point solution and one for a teit-point solution. The former should be 
adequate for all practical purposes; only where the stiffness curves are 
very irregular near the root does the ten— point solution have to be 
resorted to. The points considered by the two sets of tables are at 

= 0, 0.2, 0.4, 0.6, 0.8, and 0.9 for the shorter solution and 
8 A 

-p- = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 for the longer 
A . 

solution. The procedure to be followed for either solution is identical; 
although computing forms are presented in this paper only for the six- 
point solution, their extension to apply to the ten-point solution is 
obvious. 

Calculation of the matrices .— The first step in the aeroelastic 
analysis by means of the stiffness curves is the calculation of the 
aeroelastic matrix [A] from the physical and geometrical parameters 
of the wing. These parameters are conveniently tabulated in a form of 
the type shown in table Yl(a). The computation is then carried out 
according to the instructions of table VT(b), each step in the procedure 
being identified by the number in the upper left comer of each box. 

It must be kept in mind that many of the operations call for matrix 
multiplications where the order of the multiplicands is of importance. 

(A brief summary of matrix methods is presented in the appendix. ) The 
aeroelastic matrix is obtained as the last step (step 13) of the 
computations in this form which constitutes an evaluation of the 
A matrix given in equation (15). 
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A special case arises when e-, is zero. If e n is not zero 


matrix and the 


along the re main der of the span, its value at. some point other t.hari the 

e l / c \ 2 ^ 

root may he used as a reference value. The (—) 

[ e l r \ c r / ^ 

multiplying factors of steps 8 and 9 as well as the definition of the 
parameter a are then based on this other reference value rather 
than e l r * 3-f e^ is zero along the entire span, step 1 and steps 3 

to 8 may he omitted and steps 9 to 13 should he modified as follows : 


Step 9 

(EE) r 

q «m p 



(GJ) 

r 

®w tan A ^ 

Step 

10 

[©]• 

-[©] 

Step 

11 

As is 


Step 

12 

• 

Omit 


Step 

13 

Me r 

:0= [©] 


If influence coefficients of the proper type are available, the 
calculation of the aeroelastic matrix (A*] is carried out directly 
by means of equation (28). 

Solution for divergence dynamic pressure .— In order to dete rmin e 
the value of the parameter a or a* corresponding to divergence, the 
aeroelastic matrix QaJ or [A f ] is iterated (see appendix) as indicated 
by equation (18). Table Vll(a) may be used for this purpose. The result 
is the critical value of a or a*. The divergence dynamic pressure is 
then calculated from equation (l 6) or (27). it is to be noted that this 
pressure will be in pounds per square inch. Since the aeroelastic matrix 
is independent of the Mach numb er , except insofar as e^ varies with Mach 

number, the same critical value of a may be used to calculate the diver- 
gence dynamic pressure for an entire range of Mach numbers. If the value 
of changes, however, as it does between the subsonic and supersonic 

region, the critical value of a has to be calculated for both values 
of e-j_. 

If the value of e^ is zero along the entire span and the [a] 

matrix has been calculated according to the modified instructions, 
iteration of the matrix will give the value of the parameter d at 
divergence . From the definition of d the divergence dynamic pressure 
may then be calculated. 
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Solution for aerodynamic loading .— In order to calculate the aero- 
dynamic loading corresponding to a given. flight condition and geometric 
angle— of -attack distribution the aeroelastic matrix [A] or [A*] is 
multiplied by the value of a or a* calculated for the given flight 
condition and subtracted from the unit matrix [I]. (See equation (17).) 
The result may be entered in table Vll(b). Again it must be noted that the 
value of the aeroelastic matrix varies with the flight condition if e-^ 
varies, so that the aeroelastic matrix corresponding to the proper e^ 
value must be selected. The resulting matrix constitutes the coefficients 
of a set of simultaneous linear algebraic equations for the unknown values 

of the effective angle— of— attack distribution of the deformed wing. 
aj> in terms of the known angle— of -attack values of the rigid 

wing |ctgj. Table Vll(b) is set up for the calculation of the additional 

loading, the damping-in— roll loading, and a third arbitrary loading; as 
many loadings as desired may, of course, be calculated by this method. 

The solution of the equations may be carried out in any convenient manner. 
The form of table Vll(b) has been prepared for use in conjunction with 
Crout*a method of solving linear simultaneous equations (reference 6). 

In the case where e-^ is zero along the span, the headings at the 
top of table Vll(b) should be modified to read 



[M - aM ei=0 ] 


where [A] has been calculated according to the modified instructions 

and d has been obtained by iterating [a] q ^_q. 


The values of 


ftp) 


calculated for the additional load 


distribution (a a = l) constitute values of the ratio c-, /c^ 

/ v /, v f • *fw/ l tw 

or \ cc l/f v A CC Z/ in view of the assumptions made concerning the air 

The section loading of the flexible wing is obtained from the 


forces . 
relation 
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or in dimensionless form 


221 

c r 


m 

e l 



(32a) 


The wing lift coefficient defined by the relation 

C = — 

qS 


and the wing hending-moment coefficient defined by 


(33) 


^ (3M 

may be obtained by integrating the load distribution. These integrations 
may be performed conveniently by multiplying the cc^/c-p values by the 

first rows of the jK-J and matrices, respectively. Thus 


C Iv = 



(35) 


and 



(35a) 


_ s w c r B w 
Sj2~ 2b 



*'(36) 
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a, 


( 36 a) 
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The lateral center of pressure of the wing load 
determined from the relation 


— may then be 
s v 


y_ _ 2b Cbm^ 
^ 8 w ^^w 


( 37 ) 


The fore-and-aft location of the aerodynamic center of the wing load 
measured rearward of the leading edge of the mean aerodynamic chord as 
a fraction of the mean aerodynamic chord may be estimated from the 
relation 


(ac )^ 
C MAC 


ac + 


7 y MAC 
c MAC 


tan A ac 


( 38 ) 


where A &c is the sweep of the section-aerodynamic -center line. 

For any other geometrical angle-of-attack distributions such as 
those due to built-in twist or those due to rolling, the same section 
lift-curve slope should be used as for the structural deformations, so 

me 

that m Qi is replaced by nig and is unity in equations (32), (35), 

and (36). For the damping-in-roll distribution with a tip helix angle 
of 1 radian 



(39) 


The rolling-moment coefficient due to the wing load is defined by 


» 


2 M« 

C 7 = — ■ — £■ 

Hr qSb 


( 40 ) 


It is seen to be twice the wing bending-jHioment coefficient. 

The contribution of the wing to other stability derivatives may 
be obtained similarly by integrating the load distributions due to the 
angle-of-attack distributions caused by the motion under consideration, 
as described in reference 7; in the case of swept wings, particular care 
must be taken in selecting the proper angle-of-attack distribution nnfl 
in accounting for the lateral inclination of the lift vector. (See 
reference 4 .) 


If the aerodynamic loading or the stability derivatives are to be 
obtained for a wide variety of flight conditions, it is convenient to 
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systematize the calculations in the following manner: The aeroelastic 

matrix is computed for "both the subsonic and supersonic aerodynamic- 
center values and iterated for both cases to obtain the subsonic 
supersonic values of the divergence parameter a D . From these values 

the divergence dynamic pressure may be computed by means of equation (15) 
and plotted against Mach number, as suggested in reference 3; on the 
same plot values of the actual dynamic pressure may be plotted against 
Mach number for various altitudes of interest. Such a plot for a wing, 
the physical characteristics of which are given in figure 3, is shown in 
figure 4 . 


Since at a given Mach number the ratio a/ap is equal to the 
ratio q/q D , the range of a/a^ values of interest may’ be established 

from this plot for both the subsonic and the supersonic region. Several 
representative a/a^ values may then be chosen within the given ranges 
and the corresponding values of a computed from the previously calculated 
a p values. The aerodynamic loading is calculated for these values of a 
using the appropriate (aJ matrix and plotted in the form of 
( cc 2 ) fw/ /( cc Z . 'wi'kk "the ratio a/a^ as a parameter. From these 

curves or from the —— a values the lift coefficients may be obtained 

/ 6l 

in the form (C^f /(C]-,)^ and plotted against a/a D or q/q D ; the 
other coefficients may be obtained and plotted in a similar form. 


For any specific flight condition the value of a/a D may then 
be obtained from the plot of q and q-p against Mach number. The 

loading, lift coefficient, or other item of interest may be obtained 
from the plots which give these items in terms of the rigid-wing values. 
Once the rigid-wing values at the given Mach number are known, the 
flexible-^wing values may then be obtained immediately. 


Illustrative Example 

In order to illustrate the method described in the preceding 
sections, a typical swept wing has been analyzed. The physical and 
geometrical parameters of the wing are shown in figure 3 and the upper 
part of table VIH (which follows the form of table Vl(a)). The 
chord, eqc , and stiffness matrices have been obtained from the 
given parameters and are shown in the lower part of table VTII. 

The calculation of the aeroelastic matrix for the subsonic case 
has been carried out by means of the form of table Vl(b). All but 
three of the steps of the computation are shown in table IX numbered 
in the same order as in table Vl(b). Steps 1, 2, 6, 7, 11, and 12 
constitute matrix multiplications, which are carried out in the order 
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indicated; steps 5 and 13 constitute matrix additions or subtractions; 
steps 3 a* 1 *! 4 constitute multiplications of matrices by constants. 

The aeroelastic matrix is iterated in table X(a) (which follows 
the form, of table Yll(a)) to- yield a value of a^ = -2.208. From this 

value and a value of a^ computed in the same manner for supersonic 

speeds, the divergence dy nami c pressure has been calculated by means of 
equation (16) on the basis of estimated values of the effective lift- 
curve slope. The variation with Mach number of the divergence dynamic 
pressure, the actual dynamic pressure at sea level, and the estimated 
effective lift-curve slope is shown in figure 4. 

For a value of ■— = -0.25, such as would be obtained approximately 
a D 

at a Mach number of 1.0, the aerodynamic loading has been calculated for 
the additional— angle— of-attack case and the damping— in— roll case in 
table X(b), which follows the form of table VH(b). The values of a- 

o 

for the damping— in— roll case have been calculated from equation ( 39 ) • 

The aerodynamic loadings, in addition to those calculated for other 
q/qp values, have been plotted in figure 5 as ratios of the flexible- 
wing loadings to the rigid— wing loadings. The curves have been integrated 
to yield wing lift and rolling-moment coefficients as well as 
the aerodynamic center of the wing load, which are shown in 

table X(b) for the case of — = —0.25 and which are plotted against 

a D 

— — in figure 6. 

ID 

The wing lift coefficient is defined in such a manner that if 
the fuselage lift is known and made dimensionless by dividing by q 
and S the resulting fuselage lift coefficient may be added directly 
to the wing lift coefficient. This definition and the fact that 
figure 5(a) is plotted over the fraction of the wing-alone span Sy 

explain the fact that the area under the curve of figure 5(a) is not 1. 

The aerodynamic center as plotted in figure 6 constitutes the center 
of pressure of only the wing load. In order to obtain the airplane 
aerodynamic center, the magnitude and center of pressure of the fuselage 
load would have to be known and taken into account. 


DISCUSSION 


Both the aerodynamic and the structural assumptions made in this 
analysis are more realistic than those made in reference 3* The device 
employed in this analysis of calculating the air forces for wing sections 
parallel to the direction of flight and then transferring them to 
sections perpendicular to the elastic axis obviates the necessity of 
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replacing the actual wing with one the root and tip of which are 
perpendicular to the elastic axis for the purpose of analysis. Further- 
more, the inclusion of Falkner*s modification (see appendix) in the 
integrating matrices has the effect of rounding off the load distribu- 
tion approximately in the manner observed at subsonic speeds. At 
supersonic speeds the load distributions do not go to zero in the manner 
assumed in Falkner*s modification, but even at supersonic speeds there 
is some reduction of load at the tip, the total magnitude of which is 
not far from the reduction obtained by Falkner * s modification. 


Only one aerodynamic assumption is still made: that induction 
effects may be approximated by an over-all reduction of the strip theory 
loading (rounded off as previously described) at subcritical speeds and 
may be neglected at supersonic speeds. The effects of aerody nami c 
induction could be taken into account more accurately by using aero- 
dynamic influence— coefficient matrices instead of the effective lift— 

0 r 0 2] 

curve— elope concept and the I— 1 and matriceB used in 

L°rJ e x V°r/ 

r 

this analysis. Available methods of calculating such influence coeffi- 
cients for wings of arbitrary plan form at subsonic and supersonic speeds 
particularly those suitable for wings with large amounts of sweep, are 
either too inaccurate or too time-consuming for practical purposes, 
however. 


Although the analysis of this paper has been performed for wings 
consisting of uncambered sections, it is directly applicable as well to 
the determination of the additional loading of wings with cambered 
sections. The loading of such wings due to the section pitching moment 
at zero lift may be determined by modifying the analysis somewhat. 

The assumption of an effective root perpendicular to the elastic 
axis made in reference 3 for the purposes of calculating the structural 
response is carried over in this analysis. It is modified, however, 
to the extent that the root is no longer considered to be rigid as in 
reference 3, but flexible, both in torsion and bending. It has been 
demonstrated in reference 5 that the deflections of a swept beam may be 
estimated on that assumption, provided the root— rotation parameters are 
known. By ass uming the effective root at the intersection of the elastic 
axis with the side of the fuselage, the root bending due to bending 
moment and root twist due to torque are minimized. The bending due to 
twist and twist due to bending are the same regardless of the location of 
the effective root. 

The method of introducing the root rotations into the analysis by 
means of the [iq] matrix assures that the structural twist in planes 
parallel to the direction of flight is zero at the fuselage. From 
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figure 2 it in seen that the local values of the root rotation either 
tend to approach zero at the root or tend to cancel each other. If the 
root-rotation constants are known, the structural deformations can there- 
fore he predicted quite accurately hy the assumptions made. 

The manner in which the equations of equilibrium are solved hy 
means of the integrating matrices accounts for the true chord and 
stiffness variations. It does not necessitate replacement of the 
actual wing hy constant— chord segments with all the flexibility 
concentrated at the ends of the segments, an approach which has been 
used extensively in the work on aeroelastic problems of straight wings. 

A further refinement which obviates the necessity for making any 
structural assumptions other than that of small deflections is the use 
of measured influence coefficients in the aeroelastic analysis. Wherever 
such coefficients are available it is, of course, of advantage to use them. 

No explicit account has been taken in the analysis of the effects 
of the inertia loading on the structural deformations and hence the 
aerodynamic loading. On swept wings, in particular, their effects 
may be considerable. For the purposes of this analysis, however, the 
structural deformations due to inertia loading may be considered part 
of the geometric angle of attack and the rigid-wing geometric angle of 
attack may be modified accordingly. The deformations due to the inertia 
loading may, incidentally, be calculated conveniently by means of the E^, 
K^, and E^ matrices. 

Some of the general observations made in reference 3 concerning the 
divergence phenomenon are corroborated by the example. As expected of 
a wing with a considerable amount of sweepback, the divergence dynamic 
pressure is negative. Consequently the wing cannot diverge. The 
divergence dynamic pressure is useful as a reference value, however; 
the values of the load distribution and the stability parameters divided 
either by the corresponding rigid-wing values or by the section lift- 
curve slope depend only on the ratio of the actual to the divergence 
dynamic pressure. 

The type of plot shown in figure 4 is therefore quite, useful in 
the analysis of aeroelastic phenomena. As pointed out in reference 3 
this chart may also be used to estimate the actual divergence dynamic 
pressure where there is a possibility that the wing may diverge. It 
appears that the critical values will occur at either extremity of the 
transonic region. In the transonic region proper the lift— curve slope 
usually appears to be lower than at the extremities, so much so that the 
decrease in lift— curve slope even tends to overbalance any forward 
shift in aerodynamic center. 

As would be expected qualitatively, the effect of wing flexibility 
in the case of the example wing is to unload the wing tips owing to the 
fact that they bend up. The lift carried by the wing is therefore less 
than that carried by a rigid wing, the center of pressure being farther 
inboard and the aerodynamic center farther forward. 
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The difference between the supersonic and subsonic values of the 
loading, the lift and rolling-moment coefficients, and the aerodynamic 
center for a given value of a/a^ is due to the difference in the 

e 1 distributions. If the distributions were the same, the subsonic 
and supersonic variations of these quantities with a/a-p would coincide 
despite the difference in the e-^ values. 

Another item of possible interest is the fact that the variati ons 
of the parameters d D and a-p for the example problem are approximately 

linear (see fig. 7), as would be expected from the results of the 
analysis of reference 3. The deviations from linearity are most pronounced 
near the points for d = 0 (that is, A = 0°). They are due to the effects 
of the root rotations, in particular, the bending due to torsion and torsion 
due to bending; these effects were neglected in the approximate analysis of 
reference 3 • The points of figure 7 correspond to the example wing and 
the wings that would be obtained by rotating the example wing to the 
unswept and 37*5° sweptf orward positions in such a manner as to keep the 
e l c cos A ( EI ) 

parameters , as wel1 as the chord, stiffness, and 

moment-arm (ex) distributions constant. Points are shown for both the 
subsonic and supersonic variations as well as for the case when ep = 0 
over the entire span (a D * = 0). The difference between the subsonic and 

supersonic lines is due entirely to the difference in- the ep distribu- 
tion; if the distributions were the same, as would be the case if the 
elastic axis were at a constant fraction of the chord, the variations 
would be the same, regardless of the difference in the ep values. 

The present analysis is concerned only with wing or tail loads; 
the total loads are obtained by adding the fuselage loads (which may 
be assumed to be unaffected by flexibility) to the wing or tail loads 
obtained from the analysis. The amount of load carried by a flexible 
wing and the mann er of its distribution can consequently be estimated by 
the method presented herein if the contribution of the fuselage is known 
at low dynamic pressures, that is, for the "rigid-wing" case. 

The fuselage has a considerable effect on some of the stability 
parameters as well, although in the case of others, such as C-j , the 

‘p 

effect is negligible. Other effects that may have to be accounted for 
in calculating stability derivatives are the boundary— layer behavior 
and tip suction. The boundary— layer effect may be accounted for by 
using a section lift-curve slope corrected for boundary-layer effects 
to calculate the angle— of -attack distribution of the flexible wing at 
the flight conditions of interest and then obtaining the lift and drag 
distributions corresponding to that angle— of -attack distribution. 

Lateral tip suction may be important on low— aspect-ratio and highly 
swept wings. Since it does not affect the lift distribution, it may be 
taken into account by calculating the angle-of-attack distribution of 
the flexible wing and est ima ting the tip suction corresponding to the 
actual angle of attack at the tip. 
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lii calculating stability derivatives it is well to keep in mind 
that the method presented in this paper is based on a modified strip 
theory, unless aerodynamic influence-coefficient matrices are used. The 
calculated derivatives may therefore be somewhat in error, particularly 
if in calculating them the moment of a load distribution has to be 
determined. If there is reason to suspect that the modified strip 
theory is inadequate for calculating a given derivative, the derivative 
may be calculated for the rigid-wing case by a more refined method; the 
results calculated by the method of this paper may then be used to 
correct the accurate rigid-wing .value for the effect of structural 
flexibility. 


CONCLUDING REMARKS 


A method has been presented for calculating the aerodynamic loading, 
the divergence speed, and certain stability derivatives of wings and 
tail surfaces of arbitrary plan form and stiffness. Provisions have 
been made for using either stiffness curves and root-rotation constants 
or influence coefficients in the structural part of the analysis. Strip 
theory with overall reduction and rounding off at the tip to take 
account of aerodynamic induction have been used for the aerodynamic part 
of the analysis. Computing forms, tables of nunferical constants required 
in the analysis, and an illustrative example are included to facilitate 
calculations by means of the method. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., December 2k, 19^-8 
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APPENDIX 

SUMMARY OF MATRIX ALGEBRA PERTINENT TO THE ANALYSIS 


For the convenience of the reader unfamiliar with matrix terminology, 
a summary of matrix definitions and methods is presented in the following 
sections. For a more complete discussion of matrix methods the reader 
is referred to any text on matrices, for instance reference 8. 


Definitions 

A matrix is a rectangular array of numbers , called elements, written 
down’ in rows and columns. A column matrix consists of a single column, 
a row matrix of a single row. A square matrix has as many rows as it 
has columns. The diagonal of a square matrix from the upper left to the 
lower right is called the principal diagonal. A matrix all the elements 
of which are zero except for those on the principal diagonal is called a 
diagonal matrix. If all of these elements are unity, the matrix is 
termed the unit matrix. 


Matrix Algebra 

Addition .— Two matrices can be added or subtracted if both have the 
same number of rows and columns. The addition or subtraction is carried 
out by adding to or subtracting from each element of the first matrix 
the corresponding element of the second matrix. 

Multiplication by a constant .— A matrix is multiplied by a constant 
by multiplying each element by that constant . 

Matrix m ultiplication .— Two matrices can be multiplied by each other 
if the second has as many rows as the first has col umn s. The elements 
of the resulting matrix are obtained by multiplying the elements in the 
corresponding row of the first matrix by those of the corresponding 
column of the second matrix in the following order: The first element 

of the row is multiplied by the first element of the column, the second 
by the second, and so forth. The sum of the products obtained in this 
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manner is the value of the element of the product matrix. Schematically 
this process may he illustrated as follows: 


123456 123456 123456 



Q=aA. + bB + cC + cLD + eE + fF 

It must he emphasized that in multiplying matrices hy each other their 
order is of importance. As the two matrices under consideration are 
written the matrix at the left (the m matrix) is said to he post- 
multiplied hy the other, (the M matrix); or the M matrix may he said 
to he premultiplied hy the m matrix, in order to distinguish the manner 
in which they are multiplied. If the two matrices were written in the 
reverse order and then multiplied according to the foregoing instructions, 
that is, if the [m] matrix were postmult iplied hy the [mj matrix, the 
element of the second row and fourth column of the product matrix [M][m] 
would clearly not have the value Q in general, nor would, in general, 
any other element have the value it would have if the two matrices were 
multiplied in the order shown. Consequently it is important to observe 
the order in which the matrices are written down in the computing 
instructions. 

Matrix iteration .— The purpose of iterating a square matrix is to 
determine the col umn matrix or matrices which, if postmultiplied hy the 
given square matrix, yield the same column matrix except for a constant 
multiplier. It is the value or values of these multipliers which 
constitute the desired characteristic values of the matrix. 

The iteration is carried out hy assuming a "trial" column (the 
column shown in table VII is convenient for the purpose of this analysis) 
and premultiplying it hy the given square matrix to yield a "result" 
column. The elements of the result column including the last are 
divided hy the last element of the result column and entered as a 
second trial column. The second trial column is then premultiplied 
hy the square matrix to yield a second result column. The procedure is 
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repeated until the same value (within the desired accuracy ) is obtained 
twice in succession for the last element of the result matrix. The 
reciprocal of this value is the desired (lowest) characteristic value 
of the matrix, that is, the lowest critical value of ap, in. the 
analysis of this paper. 

Another way of estimating a first trial column in this analysis 
is to add the elements in each row of the A matrix, enter the six sums 
in the first result column, and treat them as if they had been obtained by 
multiplying the A matrix by a first trial column. 


Derivation of the Integrating Matrices 

Although familiarity with the derivation of the integrating 
matrices is not essential to the application of the method of this 
paper, an outline of the derivation is presented because of its general 
interest . 

The integrating matrices used in this paper are based on the same 
concept as Simpson* s rule — replacement of the actual function which is 
to be integrated by parabolic segments. If the function y has the 
values y n _^, y n , and y n+ i, respectively, at the equally spaced 

points Xj^, x n , and x n+ ]_, the following relations are seen to be 

true for a second— degree parabola passed through the three known points: 


y = y n + |(y n +i - y n -iX x - x n) + |(jn+i - 2 y n + y n -iX x - *n ) 2 (*i) 

f 7 “ (3 ^X-l + (3 + (5 X y n + 1 (i2 > 

(- h Ax ) 7 n-1 + (| X 7 ” + (jl (A3) 

il— 1 



r-^n+l 


y dx = 


*n 
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(x - In)y to = | to 2 )^ + (0)y n + 0 Ax^y n+1 (A5) 


J (x - xjy to = ^ ^ Jn-i + (f ^)r n -l + (S' W t" 5 > 

^TL 

where 


^ - %-l = ^n+i - % (A7) 

The different integrations over the parabolic segments may- thus 
he performed by multiplying the given y values by the multiplying 
factors indicated in equations (A2) to (a 6 ). 

Since load distributions at subsonic speeds go to zero with infinite 
slope at the tip and the ordinary second— degree parabola furnishes a 
poor approximation to such a distribution, Y. M. Falkner has suggested 
that a curve of the type 


y = A Q + A 1 (l - x ) 1 / 2 + Ag-(l - x ) 3 / 2 


(A8) 


be passed through the last three points of the load-distribution curve 
at the tip (x = l). On the basis of the approximation, relations 
equivalent to equations (Al) to (a 6 ) may be derived. The multiplying 
factors for the last two segments are then based on these equivalent 
expressions rather than those of equations (A2) to (a 6 ). 

The integrati n g factors of equations (A2) to (A 6 ) may be assembled 
directly into integrating matrices. The Ko matrix, for instance, is 

y dx. If at the upper limit x = 0.1 

0 

and if the ten-point matrix (table 111 (b)) is to be used, the factors 0.0^167, 
0.06667, and -0. OO 833 may be obtained from equation (Ak) since x^^ = 0, 

x n = 0.1, and Ax = 0.1; similarly, if for the same case the integration 
is extended to x n+ q = 0.2 as the upper limit, the integrating factors 
0*03333* 0.13333* . an( l 0.03333 will be obtained from equation (A2). 

These factors constitute the second and third rows of the matrix K 3 ; 
since the integrations are independent of the y values other than the 
first three, the other y values are multiplied by. zero in these two 
rows. In ordeivto extend the integration to x = 0.3 an integration is 



NACA TN No. 1876 


31 


again performed up to x = 0.2 and another integration, using another 
parabolic segment, is performed from x = 0.2 to x = 0.3. For the 
latter integration x^ = 0.2, = 0.3, and Ax = 0.1, .so that 

equation (A4) again yields the factors 0.014-167, O.O6667, and -O.OO833. 
The y value at x = 0.2 is therefore assigned a multiplying factor 
of O.O3333 "by the first integration and a factor of 0.014-167 by the 
second, or a total factor of 0.07500. The resulting factors are 
entered in the fourth row of the matrix. All other rows are 
obtained in a similar manner. 


The E-^ matrix is set up to perform the integration 



The 


values of the last row of the ten— point E-^ matrix (table 1(b)) are 

obtained from FaHner* s equivalent of equation (A3) for the curves 
assumed in reference 5, with = 0.8, x n = 0.9, x n+1 = 1.0, and 

Ax = 0.1. Only the m ultiplying factors for the y values at x = 0.8 
and x = 0.9 are listed, since the y value at x = 1.0 (the wing tip) 
is assumed to be zero in this analysis, so that itB multiplying factor 
is imma terial. The values of the last row but one are obtained. similarly 
from Falkner , s equivalent of equation (A2). The values of the row 


for 
x = 


JL 

S A 

0.6 


= 0.7 are obtained by using equation (A3) in the interval 
to x = 0.8 and Falkner*s equivalent of equation (A 2) in the 


interval x = 0.8 to 1.0. Similarly the row for -H- = 0.6 is 

. b A 

obtained by combining the results of equation (A2) for the interval x = 0.6 
to 0.8 with Falkner , s equivalent of equation (A2) for the 
interval x = 0.8 to 1.0. All other rows are obtained in a similar 
manner. 


The Kg matrix is set up to perform the integration 



A o 

where x is the variable of integration, and x the value of x 
the lower limit. In applying the integrating factors of equations 
through (A6) to this integration it must be realized that 


x 0 )y dx, 

at 

(A2) 



(x n - x Q ) 



(A9) 


so that the integrating factors for this integration would be obtained 

by adding (x n — x Q ) times the factors of equation (A2) or (A3) to 

the factors of equation (A5) or (A6), respectively, depending on the 

limits of the integration. The factors for the different 

segments (x = 0.8 to 1.0, 0.6 to 0.8, and so forth) are then combined 

for any given row (with its given value of x Q ) in the manner indicated for 

for the Ei matrix uo yield the E2 matrix. 
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The matrix sums up the torques outboard of a given point, 
■while the matrix gives the sum of the moments of forces applied 

outboard of a 0 given point. Neither requires any integrations in the 
sense of equations (A2) to (A 6). For the six— point method these two 
matrices are: 


*4 



• *5 


y/ ®w 

0.2 

0.4 

0.6 

0.8 

0.9 1.0 

0 

0.2 

0.4 

0.6 

0.8 

0.9 1.0 

0.2 

0 

.2 

.4 

.6 

.7 .8 

0.4 

0 

0 

.2 

.4 

.5 .6 

0.6 

0 

0 

0 

.2 

.3 .k 

0.8 

0 

0 

0 

0 

.1 .2 

0.9 

0 

0 

0 

0 

0 .1 


It will he noted that the moment arms which comprise the Ec matrix are 
fractions of s w , so that the matrix must be multiplied by the length s w 

in order to yield actual moments as stated in equation (21 ). 
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TABLE I.- VALUES OF THE INTEGRATING MATRIX 



(a) Six— Point Solution 


l/ 8 A 

0 

o2 

.1+ 

.6 

.8 

.9 

0 

0.06667 

0.26667 

0.13333 

0.26667 

0.09333 

0.15085 

.2 

-.01667 

.13333 

.15000 

.26667 

.09333 

.15085 

.4 

0 

0 

.06667 

.26667 

.09333 

.15085 

.6 

0 

0 

-.01667 

.13333 

.11000 

.15085 

.8 

0 

0 

0 

0 

.02667 

.15085 

•9 

0 

0 

0 

0 

-.01886 

.09333 


(t>) Ten— Point Solution 



0 

.1 

.2 

.3 

■1 

.5 

.6 

.7 

.8 

.9 

0 

0.03333 


0.06667 

0.13333 

0.06667 

0.13333 

0.06667 

0.13333 

0.06000 

0.15085 

.1 

-.00833 

.06667 

.07500 

.13333 

.06667 

.13333 

.06667 

•13333 

.06000 

.15085 

.2 

0 

0 

.03333 

.13333 

.06667 



.13333 

.06000 

.15085 

.3 

0 

0 . 

-.00833 

.06667 

.07500 

.13333 

.06667 

.13333 

.06000 

.15085 

.4 

0 

0 

0 

0 

.03333 

.13333 

.06667 

.13333 

.06000 

.15085 

.5 

0 

0 

0 

0 

-.00833 

.06667 

.07500 

.13333 

.06000 

.15085 

.6 

0 

0 

0 

0 

0 

0 

.03333 

.13333 

.06000 

.15085 

.7 

0 

0 

0 

0 

0 

0 

-.00833 

0 06667 

.06833 

.15085 

.8 

0 

0 

0 

0 

0 

0 

0 

0 

.02667 

.15085 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

-.01886 

.09333 
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TABLE II.- VALUES OF THE INTEGRATING MATRIX 


O] 


(a) Six-Point Solution 


n/ s A 

0 

.2 

.4 

.6 

.8 

•9 

0 

0 

0.05333 

0.05333 

0.16000 

0.07314 

0.13792 

.2 

-.00167 

.01000 

.02500 

.10667 

.05448 

.10775 

.4 

0 

0 

0 

.05333 

.03581 

.07758 

.6 

0 

0 

-.00167 

.01000 

.01548 

.04741 

.8 

0 

0 

0 

0 

-.00152 

.01724 

•9 

0 

0 

0 

0 

-.00108 

.00419 


(b) Ten— Point Solution 


Bi 

0 

99M 

.2 

.3 

.4 

\ ..5 

.6 

.7 

.8 

.9 

9 



0.013333 

0.040000 

0.026667 

0.066667 

0.040000 

0.093333 

0.046476 

0.137920 



.002500 

.006251 

.026667 



.033333 

.080000 

.040476 

.122835 

9 

0 

0 

0 

.013333 







9 

0 

0 

-.000417 

.002500 



1 

.053333 

.028476 


.4 

0 

0 

0 

0 

0 

.013333 

.013333 

.040000 

.022476 


□ 

0 

0 

*0 

0 

-.000417 

.002500 

.006251 

.026667 

.016477 

.062495 

9 

0 

0 

0 

0 

0 

0 

0 

•013333 



9 

0 

0 

0 

0 

0 

PI 


.002500. 

.004060 

— 

.032325 

.8 

0 

0 

0 

0 

0 

0 

0* 

0 

-.001523 

✓ 

.017240 

-9 

0 

0 

0 

0 

0 

0 ... 

0 



m 
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TABLE III.- VALUES OF THE INTEGRATING MATRIX 
(a) Six— Point Solution 


h/s a 

0 

.2 

a 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

.2 

.08333 

.13333 

-.01667 

0 

0 

0 

.k 

.06667 

.26667 

.06667 

0 

0 

0 

.6 

.06667 

.26667 

.15000 

.13333 

-.01667 

0 

.8 

.06667 

.26667 

.13333 

.26667 

.06667 

0 

•9 

.06667 

.26667 

.13333 

.26667 

.10833 



(t) Ten-Point Solution 



0 

.1 

.2 

.3 



.6 

.7 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

..1 

.04167 

.06667 

-.00833 

0 

0 

0 

0 

0 

0 

0 

.2 

.03333 

.13333 

.03333 

0 

0 

0 

0 

0 

0 

6 

.3 ' 

.03333 

.13333 

.O75OO 

.06667 

n 

0 

0 

0 

0 

0 


.03333 

.13333 

.06667 

.13333 

.03333 

0 

0 

0 

0 

0 

H 

.03333 

.13333 

.06667 

.13333 

.07500 

.06667 

-.00833 

0 

0 

0 

1 

.03333 

•13333 

.06667 

.13333 

.06667 



0 

0 


.7 : 

.03333 

.13333 

.06667 

.13333 

.06667 

.13333 

.07500 

.06667 

-.00833 

0 

.8 

.03333 

.13333 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.03333 

0 

.9 

.03333 

.13333 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.07500 

.06667 
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TABLE IV.- VALUES OF THE LOAD-CONVERSION MATRIX 



(a) Six— Point Solution 


n/s A 

0 

.2 

.4 

.6 

.8 

.9 

0.2 

0.08333 

0.13333 

-0.01667 

0 

0 

0 

.4 

-.01667 

.13333 

.08333 

0 

0 

0 

.6 

0 

0 

.08333 

.13333 

-.01667 

0 

.8 

0 

0 

-.01667 

•13333 

.08333 

0 

•9 

0 

0 

! 0 

0 

•04553 

.05752 

1.0 

0 

0 

0 

0 

-.01886 

.09333 




(1)) Ten— Point Solution 


l/ s A 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

0.1 

0.04166 

0.06667 

-0.00833 

0 

0 

0 

0 

0 

0 

0 

.2 

-.00833 


.04167 

0 

0 

0 

0 

0 

0 

0 

.3 

0 


.04166 

.06667 

-.00833 

0 

0 

0 

0 

0 


0 

0 

-.00833 

.06667 

.04167 

0 

0 

0 

0 

0 

D 

0 


0 

0 

.04166 

.06667 

-.00833 

0 

0 

0 

.6 

0 


0 

0 

-.00833 

.06667 

.04167 

0 

0 

0 

.7 

0 


! 0 

0 

!o~ 

0 

• 0^166 

.06667 

-.00833 

0 

.8 



0 

0 

0 

0 

-.00833 

.06667 

.04167 

0 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

.04166 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

-.00833 

.09333 
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NACA TN No. 1876 


TABLE V.- VALUES OP THE LOAD-CONVERSION MATRIX 



f 


(a) Six— Point Solution 


’l/* 'A 

0 

.2 

.4 

.6 

.8 

.9 

0.2 

0.01667 

0.16667 

0.01667 

0 

0 

0 

.4 

-.00833 

.05000 

.11667 

.05000 

-.00833 

0 

• 6 

0 

0 

.01667 

.16667 

.01667 

0 

.8 

0 

0 

-.00833 

.05000 

.08946 

.02035 

•9 

0 

0 

0 

0 

.00631 

• .08860 

1.0 

0 

0 

0 

0 

-.01077 

.04190 


(b) Ten— Point Solution 


V S A 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

0.1 

0.00833 

0.08333 

b. 00831 

V 

0 

0 

0 

p 

0 

0 

0 

.2 

-.00417 

<,02500 

.05834 

.02500 

-.00417 

0 

0 

0 

0 

0 

.3 

0 ■ 

0 

.00833 

.08333 

.00833 

0 

0 

0 

0 

0 

.4 

0 

0 

-.00417 

.02500 

.05834 

.02500 

-.00417 

0 

0 

0 

• 5 

0 

0 

0 

0 

.00833 

.08333 

.00833 

0 

0 

0 

.6 

0 

0 

0 

0 

-.00417 

.02500 

.05834 

.02500 

-. 00417 

0 

.7 

0 

0 

0 

0 

0 

0 

.00833 

.08333 

.00835 

0 

.8 

0 

0 

0 

0 

0 

0 

-.00417 

0 02500 

.06029 

.02035 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

.00631 

.08860 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

-.01077 

.04190 








































































































OF AEROELASTIC MATRIX 





























































































































































































































NACA TN No. l8j6 

:GN OP AEROELASTIC BQOATICH 


(b) Aerodynamic Loading 
































































































































































NACA TN No. 1876 


TABLE H.- COMPUTATION OP AERCKLASTIC MATRIX OF EXAMPLE WING 
(SUBSONIC CASE) 


© 

M 


l/s A 

0 

.2 

.4 

.6 

.8 

•9 

0 

0.06667 

0.21494 

0.08560 

0.13120 

0.03435 

0.04707 

.2 

— 01667 

.10746 

.09630 

.13120 

.03435 

.04707 

.4 

0 

0 

.04280 

.13120 

.03435 

.04707 

.6 

0 

0 

-.01070 

.O656O 

.04048 

.04707 

.8 

0 

0 

0 

0 

.00981 

.04707 

• 9 

0 

0 

0 

0 

-.00694 

.02912 



nn r (GJ)r 

hJ J53 T 



0 

.2 

.4 

.6 

.8 

•9 

0 

0 

0 

0 

0 

0 

0 

.2 

.08333 

.15066 

-.03501 

0 

0 

0 

.4 

.06667 

.30134 

.14001 

0 

0 

.0 

.6 

.06667 

.30134 

•31500 

.58665 

-.16070 

0 

.8 

.06667 

.30134 

.27999 

1.17335 

.64270 

0 

•9 

.06667 

.30134 

.27999 

1.17335 

1.04430 

1.04072 


© 

[ 

■a pi*] 

n/sA 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0. 

0 

O' 

0 

.2 

.08333 

.14933 

-.03201 

0 

0 

0 

.4 

.06667 

.29867 

.12801 

0 

0 

0 

.6 

.06667 

.29867 

.28800 

.49465 

-.12452 

0 

.8 

.06667 

.29867 

•25599 

.98935 

.49802 

0 

•9 

.06667 

.29867 

.25599 

.98935 

.80923 

.73070 


© 

H; 

r © 1 

i.-j 

l/sA 

0 

.2 

.4 

.6 

. .8 

•9 

0 

0 

0 

0 

0 

0 

0 

.2 

-.00613 

.01777 

.01728 

.01922 

.00503 

.00689 

.4 

-.01173 

.03724 

.04875 

.08089 

.02118 

.02902 

.6 

-.01173 

.03724 

.05082 

.17174 

.06273 

.07104 

.8 

-.01173 

.03724 

.03938 

.22118 

.10693 

.16293 

•9 

-.01173 

.03724 

.03938 

.22118 

.10254 

.23190 


® 

®- 2 a)[©] 

l/ 8 A 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

.2 

.04585 

.08216 

-.01761 

0 

0 

0 

.4 

.03668 

.16433 

.07043 

0 

0 

0 

.6 

.03668 

.16433 

.15846 

.27216 

-.06851 

0 

.8 

.03668 

.16433 

.14085 

.54434 

.27401 

0 

•9 

.03668 

• .16433 

.14085 

.54434 

.44524 

.40203 

gt-A- 0 . 55 * 


© 


*)/sA 

0 

.2 

.4 

.6 

.8 

• 9 

0 

0 

0 

0 

0 

0 

0 

.2 

-.1629 

0 

0 

0 

0 

0 

.4 

-.1629 

0 

0 

0 

0 

0 

.6 

-.1629 

0 

0 

0 

0 

0 

.8 

-.1629 

0 

0 

. 0 

0 

0 

.9 

-.1629 

0 

0 

0 

0 

0 

& (V- Q aM tenA ) = - 0 - 1629 

© 

[®M®M©] 

ti/sa 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

.2 

-.03372 

.23282 

-.05262 

0 

0 

0 

.4 

-.05955 

.46567 

.21044 

0 

0 

0 

.6 

-.05955 

.46567 

.47346 

.85881' 

-.22921 

0 

.8 

-.05955 

.46567 

.42084 

1.71769 

.91671 

0 

•9 

-.05955 

.46567 

.42084 

1.71769 

1.48954 

1.44275 


© 

Hfel 

*/*A 

0 

.2 

.4 

.6 

.8 

•9 

0 

0 

0.04826 

0.04325 

0.11456 

0.04549 

0.07917 

.2 

-.00167 

.00905 

.02028 

.07638 

.03389 

.06185 

.4 

0 

0 

0 

.03818 

.02227 

.04453 

.6 

0 

0 

-.00135 

.00716 

.00963 

.02721 

.8 

0 

0 

0 

0 

-.00095 

.00990 

•9 

0 

0 

0 

0 

-.00067 

.00241 


0 


©] [© 


*i/ 8 a 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

.2 

-.00199 

.01106 

.02444 

.08202 

.03500 

.06286 

.4 

-.00399 

.01543 

.04289 

.20681 

.09790 

.18314 

.6 

-^00399 

.01543 

.03755 

.28397 

.16542 

•33786 

.8 

-.00399 

.01543 

.03221 

.30252 

.19308 

.48336 

•9 | 

-.00399 

.01543 

.03221 

.30252 

.18680 

.52208 


© 

w-r©; 

l-[@] 

v/ B A 

0 

.2 

.4 

.6 

.8 

• 9 

0 

0 

0 

0 

0 

0 

0 

.2 

-.00414 

.00671 

-.00716 

-.06280 

-.02997 

-.05597 

.4 

-.00774 

.02181 

.OO586 

-.12592 

-.07672 

-.15412 

.6 

-.00774 

.02181 

.01327 

-.11223 

-.10269 

-.26682 

.8 

-.00774 

.02181 

.00717 

-.08134 

-.08615 

-.32043 

.9 

-.00774 

.02181 

.00717 

-.08134 

-.08426 

-.29018 
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NACA TN No. 1876 



— .Span ^ ^ 

Figure I.- Definition of geometrical parameters used in the analysi 



Rotations due to root 
deformations, <f r and T r> radians 


NACA TN No. 1876 


l 



S^= 104 in. 
w= 15 in. 

(Gj) r = 5.10 10® lb -in: 

(tlV= 3.47' 10 6 lb- in. 

T r = 43,470 in.- lb 

M r = Z 60,000 in.- lb 


Rotation 

Av. value 

Q value 

?r T 

-o.ooo£ 

0.16 

y* r M 

.0025 

.33 

Tr T 

.0010 

1.45 

* r rM 

- .0010 

- .24 



Distance along span, 17 , inches 

Figure Z.~ Rotations of q 45° swept box beam clue to root 
deflections (data from reference 5 ). 


He) to airstr 



Dynamic pressures, q and q D , pounds per square foot 

ZOO 400 600 600 1000 2000 4000 6000 8000 10000 



.4 6 .8 1.0 2 4 

Mach number, M 


Effect of Mach number on the divergence 
pressure and lift -curve slope of the example v 


Estimated effective lift -curve slope, m e , per radian 




tec,). 

(e«i) 
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NACA TN No. 1876 


*>l 

o 

o 



(a) Assumed rigid -wing loading. 



0 

- 0.25 

- 0.50 

- 0.50 


(b) Flexible -wing loading. 

Figure 5.- Load distribution of example wing- 
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delation between the parameters Q& and 6 



